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Abstract: CCD sensors do not deliver a perfect image of the light they receive. Beyond the 
well known linear image smearing due to diffusion of charges during their drift towards the pixel 
wells, non-linear effects are at play in these sensors. We now have ample evidence for both a flux- 
dependent and static image distortions, especially but not only, on deep-depleted CCDs. For large 
surveys relying on CCD sensors, these effects should now be taken into account when reducing 
data. We present here a summary of current results on sensor characterization and mitigation 
methods. 
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1. Introduction 

The first workshop on “ Precision Astronomy with Fully Depleted CCDs” was held at Brookhaven 
National Laboratory in November 2013. Most of the contributions^ discussed sensor effects affect¬ 
ing astronomical measurements using CCDs. We will here mostly summarize the findings reported 
then, and for some aspects, reported since then in the literature. A second iteration of the same 
workshop took place in December 2014 in the same place. New results where presented there, and 
we refer interested readers to the program and proceedings^. 

CCD do not map the incoming light along a perfect rectilinear and periodic grid. The depar¬ 
tures of CCD images from this ideal mapping can be broadly split into two classes: static effects that 
do not depend on the charge stored in the device, and dynamic effects that increase as charge accu¬ 
mulates in the device. A striking example of static distortions are the so-called tree rings which are 
visible as concentric variations of the sensor response to a uniform illumination. These variations 
mimic annual growth rings of trees, and qualify as static because the pattern scales with illumina¬ 
tion. We will discuss shortly the source of these patterns which were observed only recently on 
deep-depleted devices. On the contrary, dynamic distortions in CCDs vary with illumination, and 
a typical example is the variation of the size of Point Spread Function (PSF, i.e. the response to 
point sources) with intensity. Although a sizable PSF variation was only recently observed (again 
on deep-depleted CCDs), there are much earlier reports of phenomena belonging to this class: for 
example a small PSF size variation was predicted more than 40 years ago in [1]. About 10 years 
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ago, it was reported that the varianee in flat-fields does not scale with their average [2], at odds with 
Poisson statistics. This can also be attributed to dynamic distortions. 

The plan of this contribution goes as follows: We first discuss static distortions, then dynamic 
ones, and finally provide some outlook. 


2. Static distortions 


By definition static distortions are independent of the illumination level. For example a position 
dependent displacement of astronomical objects is a static distortion. A historical cause of static 
distortions of the images is due to a column or row having a physical size different from the average. 
This kind of defect might cause detectable astrometric shifts (see e.g [3], Appendix B of [4]) and 
is not specific to some kind of CCD. We do not know of such defects reported on recent CCDs (in 
particular the deep-depleted devices discussed at this workshop), perhaps because manufacturers 
have eliminated the causes. On the contrary, the defects we discuss in this section seem to be largely 
specific to deep-depleted CCDs, or at least they significantly increase with the device thickness. 


2.1 Tree rings 


The tree ring patterns of flat-field images are observed on most of the deep-depleted devices. They 
consist of concentric variations of the flat-field response to a uniform illumination. The center of 
the circles coincides with the center of the ingot from which the sensor was manufactured, as dis¬ 
cussed in [5]. On DECam, these variations amount to about 0.2% (y band) to 0.4% (g band) peak 
to peak [6], and the patterns look alike in different bands. If the flat-field image (with these rings) 
is used to actually flat-field the science images, one observes photometry offsets spatially corre¬ 
lated with the flat-field variations [6]. This suggests that these Hatfield variations do not originate 
from conversion efficiency inhomogenities but rather from pixel size variations. Indeed, applying 
to science images Hatfield variations due to pixel size inhomogeneities degrades photometric uni¬ 
formity. This hypothesis that the tree rings are due to pixel size variations is strongly supported by 
the spatial correlation of astrometry residuals and the gradient of flat-field variations [6]. 

The physical picture one can draw is the following: during the growth of the silicon boule, 
there are time-dependent variations of silicon purity, which translate into small static distortions 
of the drift field in the sensor[5, 6]. Those distortions in turn displace the image in a position- 
dependent way and the divergence of the displacement field is observed as pixel size variations 
imprinted on flat-field images. Since the tree rings h^e periods above tens of pixels, one can 
represent the displacement field as continuous. Calling 5 this displacement field, (i.e. the average 
shift in the image plane from the conversion point to the collection point), and expressing charge 
conservation, the observed Hatfield image F' reads 


F'{-t+ 5) =F{lt) 


det 


dlt 


= F{lt) ^1 — div 5 -|-o('?^)^ 


where F is the illumination pattern at the conversion point. The stronger effect at bluer wavelengths 
is naturally attributed to their average longer drift path than red wavelengths (for back-illuminated 
CCDs), and hence larger transverse shifts. This displacement field ansatz (due to drift field dis¬ 
tortions) also explains the size of the observed correlations of photometric and astrometric offsets 
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with the flat-held pattern [6]. Ineidentally, the quality of these eorrelations observed on DECam 
leaves little room (i.e. less that about one part in a thousand) for flat-field variations due to genuine 
variations of sensor effieieney, at least over spatial seales below ~ 50 pixels. 


2.2 Distorted drift fields on sensor edges 

On a CCD deviee, one might expeet distortions of the drift eleetrie field at a distanee from the 
sensor edges of the order of the sensor thiekness. For thiek deviees, this ean amount to more than 
ten pixels. Thiek deep-depleted CCDs generieally exhibit distortions of the flat-field on the sensor 
edges, either a flux defieit (a roll-off, e.g. [7]), or a flux inerease (“glowing edges” on DECam, [6]), 
depending on the CCD type. We will generieally eall roll-off the effeet in what follows, but what 
is diseussed here applies as well to brighter edge stripes. 

One ean prove that this response roll-off is due to the eleetrie field being not perpendieular 
to the sensor surfaee, again by relating astrometrie displaeements to the Hatfield ehanges. This 
is done for the e2v ESST eandidates in [7]. So, the response roll-off is another manifestation of a 
displaeement field, as tree rings, and teehniques to map tree rings will also eapture these distortions 
on the sensor edges. 

A possible handling of these defeets eonsists of just ignoring the edges. One ean note that tree 
rings eall for a general framework to handle sensor-indueed statie astrometrie distortions, whieh 
eould handle as well the response roll-off on sensor edges, in order to salvage both astrometrie and 
shape measurements. How mueh of the affeeted edge will be reeovered remains unknown today. 

The response roll-off on the edges eauses spatial eharge gradients in the image. The dynamie 
distortions we will diseuss in next seetion tend to smear eontrasts, more speeifieally to eause larger 
smearing to larger eontrasts. As a eonsequenee the shape of the roll-off does not exaetly seale 
with illumination, as was observed on two ESST sensor eandidates [8]. For seienee observations, 
this means that the handling of the edge roll-off should in prineiple depend on the sky baekground 
eolleeted in the image. However, if one implements the eorreetion of dynamie distortions at the 
pixel level (as suggested in [9, 10]), the dependenee on the sky baekground of the response roll-off 
is aeeounted for in this eorreetion. The bulk of the effeet still remains to be eorreeted for, using 
some mapping of the displaeement field. 


2.3 Static distortions and flat fielding 

When redueing CCD data, it is fairly eommon to divide every image by some average of flat-field 
frames in order to eompensate for effieieney variations. When the spatial variations of flat-field 
are dominated by pixel size variations, this is likely to degrade the photometrie uniformity [11]. 
Assuming the eontribution of pixel size variations to the flat-field ean be evaluated, one might be 
tempted to take those out of the flat-field. This would leave spatial variations of the sky baekground 
in seienee images whieh would make the sky subtraetion inaeeurate if not unpraetieal. So, over 
spatial seales smaller than the ones used for sky baekground evaluation, the pixel size variations 
have to be left into the flat-field frame and their effeet on photometry should be eorreeted for in 
image eatalogs. 
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3. Dynamic distortions 


Dynamic distortions are characterized by departures from a linear response: the relation between 
the image and the illumination (slightly) depends on the illumination level, even well below the 
saturation level of the sensor, and independently of possible non-linearity of the electronic chain. 

3.1 Fat PSF or the brighter-fatter effect 

The growth of the PSF width with flux (e.g. [9]) is a striking example of dynamic response dis¬ 
tortion. The increase in size reaches ~ 2% from zero flux to saturation on one of LSST candidate 
sensors (the CCD 250 from e2v, 100 jJ. thick, 10/tm pitch). The size increase is smaller (~ 1%) 
on the DECam sensors (LBNL CCDs, 250 /tm thickness, 13.5 /tm pitch). We do not know of any 
failed attempt to observe this effect, even on thinned CCDs ([4]), where is at only a few per mil 
level. We do not know either of any reported evidence of chromaticity of the effect. All reports are 
compatible with a slightly steeper increase (by ~ 20 %) of the PSF size along the columns than 
along the rows. This phenomenon qualifies as non-linear because the response does not exactly 
scale with illumination. 

A relative increase of 1% of the PSF size should be properly addressed for at least two sci¬ 
entific applications of astronomical imaging : precision PSF photometry of faint vs bright sources 
(necessary for example for measuring luminosity distances to supernovae), and cosmic shear mea¬ 
surements via galaxy ellipticities (which involve comparing the second moments of faint galaxies 
and brighter stars). Qualitatively, the brighter-fatter effect is a smearing of the image increasing 
with illumination, and the next section deals with another behavior that could be described using 
the same expression. 

3.2 Departure from Poisson fluctuations in uniform exposures 

The non-linearity of the Photon Transfer Curve (PTC) (e.g. [2]) is some other unexpected behavior 
of CCDs. Namely one observes that the spatial variance of uniform illuminations does not scale 
with their average, but slightly falls off, at odds with expectations from Poisson statistics. Unsur¬ 
prisingly, this “variance deficit” comes with (mostly) positive correlations between nearby pixels 
which decay with spatial separation. These correlations also scale with the average illumination 
[2, 9, 12]. 

The fact that the variance does not scale exactly with the average qualifies indeed as non¬ 
linearity, as the following sketch shows. Fet us split a uniform long illumination into smaller time 
chunks. If the device is perfectly linear, the various chunks are independent and both their averages 
and variances just add up when assembling the long exposure. So, a non-linear PTC indicates a 
response non-linearity, and one can suspect that because of causality, when integrating an image, 
the late stages of charge accumulation are influenced by the outcome of earlier stages. 

So, the bending of the PTC is non-linear, and corresponds, as the brighter-fatter effect, to a 
smearing of the response, increasing with the exposure level. Since we observe that the statistical 
correlations, up to sensor saturation, grow as the average jx of the Hatfield, we get that the covari¬ 
ance of pixels separated by i rows and j lines vary like Covij jlV, where V is the variance of the 
Hatfield. In a charge-conserving process, this evolution of covariances comes with a compensating 
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Figure 1. Electrostatic computation of field lines in a 100 /rm thick, 10 /rm pitch CCD. The clock lines are 
perpendicular to the picture plane, the picture only displays the 20 /rm part that contains the clock stripes, 
and the buried channel has not been implemented in the setup . One set of field lines applies to the empty 
device, and the other applies to the device with 50 000 ke in one of the potential wells (red spot). In both 
instances, the thick field lines are the effective pixel boundaries. Note that adding this charge actually shrinks 
the pixel in which it resides, and also to a lesser extend its neighbors. Note that the effect of the distortions 
is concentrated in the last microns of the drift and hence essentially do not depend on the conversion point 
of the photons. The dynamic distortions are then mostly achromatic, as observed. This figure is borrowed 
from [13]. 

evolution of the variance (which is hence not exactly Poissonian), because the integral of the cor¬ 
relation function is conserved by charge redistribution (see e.g. the appendices of [2, 13]). In both 
[2, 13], it is also checked that summing the correlations actually restores the linearity of the PTC. 

With Covij oc IJ.V, i.e. Cov,/ oc to first order of perturbations, one expects that the variance 
of uniform exposures is a quadratic function of their average, as observed (e.g. [2, 9, 13]). Note 
that the quadratic correction to Poisson statistics is not dominant: the variance deficit is typically 
of the order of 10 % at saturation. The correlations between nearest neighbors are in the percent 
range at saturation and rapidly decay to typically 10^^ at a few pixel separation. 

3.3 Relating correlations and the brighter-fatter effect 

In [9] it is proposed that correlations in uniform exposures and the brighter-fatter effect share the 
same origin, namely the distortions of the electric field sourced by fhe charges stored in fhe pofen- 
fial wells of fhe sensor, as illusfrafed by an elecfrosfafic compufafion in Fig. [I]. This assumption 
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is supported by crude electrostatic calculations, which show that adjusting the size of field dis¬ 
tortions induced by some arbitrary charge, one can reproduce (to significantly better than a factor 
of 2) both the flatfield correlations and how rapidly spot widths increase with flux. Esfablishing 
how sfored charges disforf fhe field lines in fhe vicinity of fhe charge buckefs requires a precise 
knowledge of fhe geomefry and nafure of fhe clock sfripes and channel slops. This is usually nol 
provided fogefher wifh sensors, and some of fhis informafion is even regarded as propriefary by 
sensor vendors. If is fhus fempfing fo derive how sfored charges affecl drifl palhs from measured 
correlalions in uniform exposures, in order lo infer fhe brighler-faller effecf. Note lhal fhis simple 
physical model of disforfions also explains why fhe observed effecls are mosfly achromatic: sfored 
charges affecl fhe eleclric field mosfly in fhe vicinity of fhe pixel wells (see Fig. |T]), and hence 
affecl all wavelenglhs equally. We mighl however expecl some manifeslalions of chromalicify of 
correlalions al large dislances. The reduction of fhe drifl eleclric field by charges stored in fhe 
polenlial wells increases laleral diffusion (see e.g. [14] and references Iherein), and fhis increase is 
generally a small conlribulion lo fhe brighler-faller effecf [13], bul if is chromalic as well. 

If fhe brighler-faller effecl manifesls ilself as a moslly linear increase of Ihe PSF as a function 
of flux (or peak flux), one cannol exclude lhal Ihe effecl is more lhan jusl a scaling of size. Iden¬ 
tifying fhe physical mechanism lhal causes Ihis apparenl increase is Ihe way to safety model Ihe 
average evolution of Ihe shape of poinl-like objecls as a function of Iheir flux. Once Ihis is done, 
brighl slars can safety be used to e.g. unfold Ihe finite PSF size from measuremenls of fainl galaxy 
shapes, as required by cosmic shear measuremenls. 

Ref. [9] proposes a lirsl order empirical paramelrizalion of pixel boundary shifts induced by 
stored charges. This simple model describes Ihe correlations in uniform exposures linearly rising 
wilh Ihe average conlenl, Ihe quadratic behavior of Ihe PTC^, and Ihe essentially linear rise of Ihe 
PSF size wilh peak flux. The critical lesl of Ihe model (in facl of any model relating correlations 
and Ihe brighler-faller effecl) consisls in exlracling Ihe model parameters from flal-fields, and Ihen 
quanlilalively compare Ihe predicted PSF widlh evolution wilh flux wilh Ihe measured one. 

The model proposed in [9] has aboul Iwice more parameters lhan available correlations to 
measure. Ref. [9], [13], and [10] follow similar slralegies in order to overcome Ihis limitation: 
they assume that the electrostatic forces sourced by stored charges decay smoothly with distance. 
The details of how smoothness is imposed do not seem to alter significantly the result. These three 
works provide evidence that the brighter-fatter slope can be predicted from correlations in uniform 
exposures. Ref. [13] find a 5% mismatch compatible with the statistics in use. 

4. Further work and outlook 

We seem to have now physical explanations for both static and dynamic distortions: the former are 
due to the drift electric field (in the bulk of the device) not being orthogonal to the surface, either 
because of inhomogeneities of the bulk properties, or to the electric field being distorted on the 
edges. The dynamic distortions are due to the electric field being affected by the charges stored 
in the potential wells. Because the dynamic distortions mostly happen near the charge buckets, 
and the static ones along the whole drift path, static distortions are chromatic while the dynamic 

^The linear increase of correlations in uniform exposures implies a quadratic variation of the PTC, as soon as charge 
is conserved. 
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ones are mostly achromatic. Note that although the edge roll-off is both static and dynamic, the 
separation we are proposing remains pertinent. 

Static distortions are important to map out because they affect both flux and shape measure¬ 
ments. A detailed study of the static distortions in DECam is presented in [6] and proposes a simple 
recipe to handle the tree rings: it determines the displacement field from flat-field inhomogenities, 
assuming that these perturbations are entirely due to displacements, and that displacements follow 
a radial pattern. This last conditions allows one to turn a scalar field (the Hatfield inhomogenities) 
into a 2-d vector displacement field in the image plane. Whether this simple method can be applied 
in general is questionable, and a practical method to separate quantum efficiency variations from 
a displacement field is still to be demonstrated. Since the 2-d displacements are due to transverse 
drift fields, one can argue that this 2-d displacement field has a null rotational, which turns it back 
into a scalar field. However, one would probably keep this property for a sanity check rather than 
reducing the solution space to fields deriving from a potential. Incidentally, cosmic shear fields 
(which are as well irrotational) are usually not evaluated under this constraint (see e.g. [15]), but it 
is rather used as a check. 

In order to separate pixel size effects for genuine sensor efficiency variations, we do not know 
of a simple method so far. We can think of two avenues. In the laboratory, one can illuminate 
devices with spatially varying patterns, and with enough rotations and shifts of the same pattern, 
one can solve for the displacement field, sensor efficiency, and, if needed, the illumination pattern 
itself. Several enterprises along these lines are being developed but the proof of concept is still 
to be delivered. On the sky, one can measure the displacement field from astrometry (e.g. [6]), 
evaluate the effect on flat-fields and attribute what is left in flat fields to sensor efficiency variations. 
Reasonably, projects now assembling their hardware (like LSST) should follow both avenues. 

Reference [10] constitutes arguably the first large scale attempt to infer the brighter-fatter effect 
from pixel correlations, and the results are encouraging. One might then regard electrostatic mod¬ 
elling as a useless complication, especially since this work detects correlations slopes that change 
from sensor to sensor, which are difficult to reduce to variation between sensor batches. In our 
opinion, it is still pertinent to develop full realistic electrostatic models, because we cannot exclude 
that the sensor parameters required for these simulations might be constrained from measurements 
of correlations. Furthermore, since we expect that some level of chromatism of correlations is 
eventually detected, this might further constrain the simulated electrostatic configurations. 
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